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any  given  frequency  there  always  exist  some  set  of  structural  wave  numbers 

(Jtn,kn)  that  represent  freely  propagating  waves.  This  result  is  in  contrast 
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to  previous  one-dimensional  analyses  that  result  in  specific  tests  of 
propagating  and  nonpropagating  bands. 


The  nature  of  the  response  of  the  structure  as  a  function  of  frequency  is 
provided  in  terms  of  the  various  inertial  and  elastic  properties  of  the 
structure,  and  an  analogy  is  drawn  with  the  one-dimensional  result  in  terms 
of  propagating  wave  number  sets  being  confronted  with  an  inertial  or 
compliant  impedance .  The  analogy  results  in  the  propagating  pass  bands 
flipping  about  the  one-dimensional  dlamped-clamped  resonant  frequencies.  _ 

Numerical  results  are  presented  comparing  the  on-axis  radiated  pressure 
associated  with  the  one-dimensional  and  two-dimensional  grid  structures. 

For  the  drive  point  at  a  beam  location  the  difference  in  on-axis  radiated 
pressure  is  due  to  the  larger  inertial  impedance  of  the  two-dimensional 
grid  structure. 
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ABSTRACT 

A  general  solution  has  been  developed  for  the  transform 

response  and  the  radiated  pressure  associated  with  a  fluid  loaded, 

orthogonal  grid  supported  plate  driven  by  a  harmonic  point  force. 

For  the  case  of  light  or  no  fluid  loading  the  solution  is  in  a 

closed  form.  A  particular  interest  in  developing  this  solution 

was  to  determine  the  structural  pass-band  characteristics  of  such 

a  structure  in  terms  of  the  inertial  and  elastic  properties  of  the 

beams  composing  the  grid  as  well  as  the  plate  itself.  It  is  shown 

that  at  any  given  frequency  there  always  exist  some  set  of  structural 

wave  numbers  (k  ,k  )  that  represent  freely  propagating  waves.  This 
x  y 

result  is  in  contrast  to  previous  one-dimensional  analyses  that 
result  in  specific  sets  of  propagating  and  nonpropagating  bands. 

The  nature  of  the  response  of  the  structure  as  a  function  of 
frequency  is  provided  in  terms  of  the  various  inertial  and  elastic 
properties  of  the  structure,  and  an  analogy  is  drawn  with  the  one¬ 
dimensional  result  in  terms  of  propagating  wave  number  sets  being 
confronted  with  an  inertial  or  compliant  impedance.  The  analogy 
results  in  the  propagating  pass  bands  flipping  about  the  one¬ 
dimensional  clamped-clamped  resonant  frequencies. 

Numerical  results  are  presented  comparing  the  on-axis  radiated 
pressure  associated  with  the  one-dimensional  and  two-dimensional 
grid  structures.  For  the  drive  point  at  a  beam  location  the 
difference  in  on-axis  radiated  pressure  is  due  to  the  larger 
inertial  impedance  of  the  two-dimensional  grid  structure. 
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I .  INTRODUCTION 

One  of  the  basic  building  blocks  of  hull  or  cabin-like  structures  is 
the  plating-stringer  combination.  Typically,  machinery  induced  vibratory 
forces  are  applied  to  a  stringer  through  a  foundation  structure.  The 
vibration  field  then  propagates  along  the  hull  via  the  hull  plating  as 
well  as  the  stringers. 

1  2 

Cremer  and  Heckl  and  Mead,  have  analyzed  the  one-dimensional  case, 
i.e.,  wave  propagation  in  periodically  supported,  undamped  beams.  They 
showed  that  within  discrete  frequency  bands  waves  can  propagate  freely 
without  decay.  Outside  these  "propagation  bands,"  waves  decay  with 
distance.  The  particular  problem  considered  by  Mead  assumed  that  the 
periodic  supports  offer  infinite  translational  impedances  and  spring-like 
rotational  constraints.  Under  the  assumption  that  the  excitation  is 
localized  to  within  a  single  bay,  he  found  that  wave  motion  is  produced  in 
the  adjacent  bays  and  that  these  waves  are  freely  propagating  within 
discrete  frequency  bands.  The  purpose  of  this  report  is  to  analyze  the 
cut-off  frequencies  and  pass-bands  associated  with  the  two-dimensional 
problem,  namely  a  point  driven  orthogonal  grid  structure,  and  to  interpret 
those  results  in  terms  of  the  propagation  of  structure-borne  noise  and  the 
radiation  of  sound. 


U 
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II.  ANALYSIS 

A.  General 

The  structure  to  be  undertaken  in  this  analysis  is  that  of  an 
infinite  plate  periodically  supported  by  a  square  orthogonal  grid  and 
submerged  in  a  semi-infinite  acoustic  medium.  The  problem  to  be  analyzed 
is  the  response  of  this  structure  to  a  single  harmonic  point  force  applied 
to  the  plate  at  a  grid  intersection  point.  Specifically ,  we  will  examine 
"cut-off"  phenomena  in  terms  of  the  various  elastic  properties  of  the 
structure.  The  subsidiary  result  for  the  case  of  a  one  dimensional 
periodic  array  of  beams  supporting  the  plate  is  also  provided. 

The  geometry  utilized  in  the  analysis  is  depicted  in  Fig.  Il-la  and  b. 

B.  Equations  of  Motion 

In  the  following  development  we  have  assumed  a  harmonic  time 
dependence  of  the  form,  e  *u,t.  Derivatives  with  respect  to  time  in  the 
plate  and  beam  equations  have  thus  been  taken,  and  the  time  dependence 
factored  out. 

1.  Thin  Plate  Equations  -  Assuming  that  the  force  is  applied  at 
the  origin  which  is  also  a  grid  intersection  point,  then  the  equation  of 
motion  is, 


[DV4-m  ui2]w  (x,y)  =  F  S(x)  6(y)  +  l  Q(x,y)  S(x-nL) 
p  p  o  n*-oo 


+  1  Q(x,y) 6(y-mL)  +  p(x,y) 

m=-°° 


where 


D  -  E  h3/12(l-v2) 
P  P  P 


Ep  »  Plate  modulus  of  elasticity 


h  -  Plate  thickness 
P 


v  *  Poisson's  Ratio 
P 


(II-l) 
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m  •  p  h 

P  P  P 

pp  ■  Density  of  plate 

*4  2 . 

k_  -  m  w  /D 
f  P 

Fq  »  Applied  Force 

Q{x,y)  »  Beam  reaction  force 

p(x,y)  »  Fluid  pressure  reaction  on  plate 

Upon  applying  the  double  Fourier  transform  to  Eq.  II-l,  we  obtain 


■  LF 


o*J-.  »[("*  -  T1)  •  s] 
«[v  (*y  -  *r)]  ♦ 


( 1 1 — 2 ) 


+  Z 
m=-<» 


k  ) 

y 


where  we  have  utilized  the  following  transform  relations, 

-ik  x  -ik  y 


CD  CD 


f (k  ,k  ) 
x  y 


II  f(x,y)e  X  e  y  dxdy 


—  oo  —CO 


OO  CO 


f(x,y) 


Iff 


„  +ik  x  +ik  y 

f(k  ,k  )  e  X  e  y  dk  kd 

x  Y  x  y 


— oo  — oo 


/ 


-ik  x  ce 

f(x)5(x-nL)e  X  dx  *»  —  E  f 

L  n—* 


(*x  -  ¥) 


Now  for  a  plate  of  infinite  extent  we  can  also  write  the  transformed  surface 
pressure  in  terms  of  the  transformed  plate  displacement. 3 
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where 

p  -  Density  of  surrounding  fluid 
k  -  Acoustic  wave  number  in  fluid  (u>/c) 


2.  Beam  (Grid)  Equations  -  We  can  write  two  equations  for  the 
grid  network  of  orthogonal  (Euler-Bernoulli)  beams, 


3  2 

fEBIB  — 4  ~  V  lwD(x,y)  S<x-nL)  ■  -  Q(x,y) 5(x-nL) 

ay  P 


(II-4) 


8  2 

IEb:b  — 4  ”  mbu,',w n(x'y)  5<y-mL>  “  -  Q(x,y)  <5(y-mL) 
3x 


where 


Eg  -  Beam  modulus  of  elasticity 


Ifi  ”  Beam  moment  of  inertia 


Mg  *  PgAg,  mass  P«r  unit  length  of  beam 


( II-5) 


V  /EBrB 


Applying  the  double  Fourier  transform  to  the  two  beam  equations  results  in, 

i»W,4p[(k«  -  ¥)  •  \]  -  -  «[(kx  -  ¥ )  •  s] 


Vb- 


4  4  ^ 

El  (x  -k|!) w 
B  B  x  B  p 


Mky-¥)]--«[v( 


ky  -  ¥)] 


(II-6) 


(II-7) 


C.  The  Coupled  Solution 

In  developing  the  solution  to  the  coupled  equations  of  motion 
Eqs.  II-2,  II-6,  and  II-7)  we  have  defined  the  following  admittance 
functions. 
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I  w  (k*,k  )  *  F  I  Y  (k*,k  )  +  I  Q(kn,k  )  Z  Y  (k*,k  ) 
^  P  x  y  o^pxy  n  xy^pxy 


+  Z  Y  (k*,k  )  Z  Q(k\kn) 
jpxyn  xy 


(11-13) 


Similarily  using  Eq.  11-10  again  but  now  setting  k^  =  k^  and  summing  over 
all  i,  we  obtain. 


Z  w  (k  ,kS  =  F  Z  Y  (k  ,k*)  +  Z  Y  (k  ,kS  Z  Q(kn,kS 
4pxy  °  ^  P  x  y  p  x  y  n  *  x'  y 


(11-14) 


♦  s  °(vk">  \  ykX> 

In  both  of  these  manipulations  we  have  utilized  the  following  identity 

Z  A(n)  l  C (n+m)  =  Z  C(m)  Z  A(n) 
n  m  m  n 

Now  summing  over  Eq.  11-11,  and  11-12  we  have 


Z  w  (kn,k  )  =  -Y*  Z  Q(kn,k  ) 
n  p  x  y  B  n  x  y 


(11-15) 


Z  w  (k  ,kn)  =  -Y*  Z  Q(k  ,kn) 
npxy  Bn  x  y 


(11-16) 


Substituting  Eq.  11-15  into  11-13,  and  11-16  into  11-14  yields  the  following 
set  of  equations, 


£  Q(k",k  )  =  Fo  E  Y  (k\k  )  +  Z  Q(kn,k  )  z  Y  (k4,k  ) 
Bn  x  y  o  p  x  y  n  xy  pxy 

+  Z  Y  (k*,k  )  Z  Q(k£,kn) 
^Pxyn*x  y 


(II-18) 


-Y*  Z  Q(k  ,k*>  =  F  Z  Y  (k  ,k^)  +  Z  Y  Ck  ,k4)  Z  Q(k"  k4> 


x  y 


°  n  P  x  y 


p  x  y  n  ~  x  y 


+  Z  Q(k  ,kn)  Z  Y  (k  ,k*) 
n  x  y  pxy 


-6- 


(11-19) 


*  't'V - r 
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To  solve  for  the  line  reaction  transform  in  terms  of  the  known  admittance 
functions  we  must  make  use  of  the  symmetry  inherent  in  the  problem.  Thus, 
by  symmetry  we  know. 


S  »"vV  ■  5  “‘’VV 


(11-20) 


Multiplying  both  sides  by  Y^tk^k^)  '  rePlacing  with  k^  and  summing  over 
all  m  results  in  the  following  identity. 


E  Y  (k  ,k  )  T.  Q(k  ,k  )  =  I  2(k  ,k  )  EY  (k  ,k  | 
fflpxyn  xy  n  xympxy 


(II-21) 


In  a  parallel  development  we  can  show  that 


E  Y  (k  ,k  )  £  Q(k  ,k  )  =  E  Q(k  /k  )  E  Y  (k  ,k  ) 
npxym  xy  m  xynpxy 


(11-22) 


Substituting  these  two  results  in  Eq.  11-18,  and  11-19  respectively,  yields 
the  following  pair  of  equations. 


-  E  Q(kn,k  ) [Y^  +  S  ]  =  F  S  +  S  E  Q(k  ,km) 
n  x  y  B  x  ox  ym  xy 


(11-23) 


- E  Q(k  ,km)  [Y*  +  S  ] 
m  x  y  B  y 


=  F  S  +  S  E  Q(k  ,k  ) 
o  y  x  n  x  y 


(11-24) 


where 


S  =  E  Y  (k  ,k  ) 
x  n  p  x  y 


(11-25) 


S  =  E  Y  (k  ,k  ) 
y  m  p  x  y 


(11-26) 


Solving  Eqs.  11-23  and  11-24  simultaneously  for  the  transform  grid  reaction 
sums  and  substituting  into  Eq.  11-10  we  reach  the  final  result  for  the 
transform  of  the  displacement  response. 


w  (k  ,k  )  =  Y  (k  ,k  ) F 
p  x  y  p  x  y  o 


Y*S  +Y*S 
B  X  B  \ 


(11-27) 


I 
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I 


D. 


The  Plate  Admittance  Sums, 


S  and  S 
_* _ 1 


The  result  provided  by  Eq.  11-27  is  limited  to  numerical  computation 
by  the  infinite  sums  of  the  plate  admittance  functions.  Examining  these  sums 
we  have  from  Eq.  11-25, 


S 

x 


Z 

n“-° 


Y  (k  ,k  ) 
p  n  y 


S 

y 


? 

m=-“ 


Y  (k  ,km) 
p  x  y 


where. 


(Y  (kn,k  )]_1 
p  x  y 


,2  212 

4  ) 

+  Kv 

"  kf  I 

1  yJ 

) 

J 

iaipL 

( 

<  [*a  -  (k*  - 

¥)2- 

^  ] 1/2  * 

(11-28) 


We  can  see  that  with  fluid  loading  included  this  sum  will  be 


2  2  2 
values  of,  k  <  k  +  k  . 

x  y 


2  2  2 

For  values  of  k  >  k  +  k  the  sum 
-  x  y 


both  real  and  imaginary  components. 


wholly  real  for 
will  result  in 


With  the  fluid  loading  the  sums  must  be  evaluated  numerically.  There  is 
no  problem  with  convergence  but  as  higher  frequencies  are  considered  more 
terms  will  be  required.  In  general  the  number  of  terms  needed  is  the  order 
of  2kfL. 

Neglecting  the  effect  of  fluid  loading,  however,  the  sums  are  in  the 
form  of  a  quartic  in  n. 


S 

x 


(i £)(£ i 


V"'VY 


(11-29) 


where. 


p  (n,f  ,f  ) 
x  x  y 


.2  .2, 2  ,4 

Kfx-„>  ♦  y  -  fp 


k  L 
x 

2tt 


(II-30a) 


(I I- 30b) 


at* 
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k  L 


(II- 30c) 


Now  the  quartic  in  the  denominator  of  the  sum  has  four  roots 


Px(n,fx,fy)  »  (n-Pj^)  (n-p2)  (n-p3)  (n-p^) 


( II- 30d) 


(11-31) 


where, 


=  f 

v2  x 


1/2 


1/2 


f  +  i(f2  +  fM 
x  V  P  Y  / 


£  -  1  (£p  +  fv) 


1/2 


1/2 


*4 

Therefore  using  the  residue  theorem  result 

oo  . 

z  p(n)  =  -  (sum  of  the  residues  of  TTp(z)cot"z  at  the  poles  of  p(z)  ) 


we  can  write  the  result  for  Sx  without  fluid  loading. 


S 

x 


\LD/  \2rr /  \4f2^ 


cot7r[fx->-(f^-fV/2)-cotTr[fx-(fp-fy/2) 


(f2-f2)1/2 

p  y 


(11-32) 


icotw  lfx+i  ( f 2+f y)  i//2)-icotir  [ f x~ i  ( fp+fy)  1 

«y//2 


The  result  for  Sy  is  directly  obtained  by  simply  switching  fx  and  fy  in 
Eq.  11-32. 


-9- 
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E.  The  Transform  Displacement  Response 

Using  Eq.  11-29  with  the  results  just  produced,  we  can  write  out 
the  transform  of  the  displacement  response  in  a  closed  form  with  no  infinite 
sums.  The  most  general  expression  does  not  appear  to  be  simply  amenable  to 
inverse  transformation,  however,  we  can  obtain  the  information  on  passbands 
from  the  behavior  of  the  transform  itself.  To  do  this  it  will  be  simple  and 
more  enlightening  to  examine  our  general  expression  in  some  of  the  limiting 
regions  of  wavenumber  space  where  the  algebraic  form  of  the  result  is  not 
quite  so  ponderous.  These  results  are  undertaken  in  the  balance  of  this 
section. 


1.  Case  #1,  k,  >>  k  ,k  -  In  this  case  the  plate  admittance  sums 
f  x  y 


are  independent  of  k^  and  k^. 


S  =  S 
x  y 


-kf  r  kfL  kfL- 

- 9  cot  +  coth  ~T~ 

4m  a)  _ 

n  L 


(IT-33) 


y  (k  ,k  )  =  - 
P  x  y 


Lm  ui 
P 


(11-34) 


Assuming  the  beam  wavenumber,  k  ,  to  be  of  arbitrary  magnitude,*  we  can  write 

B 


(11-35) 


w 


*  For  typical  beam/plate  combinations  the  ratio  of  k_/k.  would  be  expected  to 

b  r 

be  less  than  unity. 
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2.  Case  #2,  kf  »  k  ,k  and  k  >>  k  ,k  -  For  this  case  we  can 

_ £ _ x  y _ s _ x  y 

write  the  answer  directly  from  our  previous  result,  or 


w  (k  ,k  ) 
p  x  y 


V  fi  + Y 


.(¥)(-¥•-*¥) 


(11-36) 


Note  that  Eq.  11-36  is  also  the  solution  for  the  very  important  case  of 

k  »k  =  0.  This  solution,  related  to  the  on-axis  radiated  sound  pressure, 
x  y 

is  discussed  in  greater  detail  in  Section  III. 

3.  Case  #3,  k  »k  >>  k  ,k  ,  and  k  L,k  L  >>  1  -  In  this  instance 

_ x  y _ r  s _ x  y _ 

we  have  used  the  results, 

2 

Lim  esc  (x±iy)  ■  0 
x,y-x» 

Lim  cot(x±iy)  -  +  i 

x,y^» 

to  arrive  at  greatly  simpl*fied  expressions  for  the  plate  admittance  functions, 


x  ,  3 

40  k 

y 


i 

4D  k  3 
x 


and  our  resultant  displacement  response  is, 


w  (k  ,k  ) 
P  x  y 


2  2  2  2 
mu  (k  +k  ) 

P  x  y  1  + 


t(c) 


L(  k  +  k  ) 

x  y 


(11-37) 


4.  Case  #4,  k  ■  k  «  k  -  This  particular  case,  of  course,  is 

_ _  _  Y  i 


limited  to  a  single  point  instead  of  a  region  in  wavenumber  space.  There  are 
two  solutions  with  simplified  results. 
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a.  k^L  >>  1 


w  (k  ,k  ) 
p  x  y 


3m  a) 
P 


Y  /kr 

l\  r  0 

1  +  — f-f 
2/2  '  2 

-)  [*  V  ~2  (“f“) 

-fc)  -iJ 

(11-38) 


b. 


w  (k  ,k  )  = 
p  x  y 


«  1 


3m  u 
P 


_ l 


(11-39) 


F.  One  Dimensional  Beam  Array 

A  subsidiary  result  of  the  preceding  analysis  for  the  grid  supported 
plate  is  the  solution  for  the  point  driven  plate  supported  by  a  one  dimensional 
periodic  array  of  beams  (Fig.  Il-lb) . 

As  the  result  follows  directly  from  the  grid  analysis  we  will  simply 
outline  the  steps  beginning  with  the  transformed  equation  of  motion. 


Plate 


y  ^(k  ,k  ) w  (k  ,k  ) 
P  x  y  p  x  y 


F  +  E  Q(kn,k  ) 
o  n  *  x  y 


(11-40) 


Beam 


y”L(k  )w  (kn,k  )  =  -  Q(k*,k  ) 
B  ypxy  ny 


(11-41) 


c 


where  y^1  and  y^1  are  given  by  Eqs.  II-8  and  11-9.  Using  the  same  approach  as 
followed  in  developing  Eqs.  11-13  through  11-16,  we  can  reach  the  following 


result 


-F  S 

l  Q(kn,k  )  =  * 

n  x  y  yy+s 

B  x 


(11-42) 


1 
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and  the  transform  displacement  is 


w  (k  ,k  ) 
P  x  y 


y  (k  ,k  ) F 
P  x  y  o 


1  + 


yn(k  ) 
B  y 


(11-43) 


where  is  still  provided  by  the  result  shown  in  Eq.  11-32.  As  with  the 
grid  solution  we  can  develop  the  following  limiting  results. 


1.  Case  #1,  k.  >>  k  ,k 
_ f _ x  y 


w  (k  #k  ) 
p  x  y 


m  w2  )  Y 
P  (.!--?■ 


kfL 


'k 

JL 

kB> 


-  1 


>  (11-44) 


kfL  k  L 

cot  — —  +  coth  - 

2  2 


2.  Case  #2,  k.  >>  k  ,k  and  k  >>  k 
_ f  x  y _ B  y 


w  (k  ,k  ) 
P  x  y 


Zo_  \ _ i _ l 


(11-45) 


This,  of  course,  is  also  the  result  for  k  =  k  *=  0. 

x  y 


3.  Case  #3,  k  ,k  >>  k.,k  and  k  L,k  L  >>  1 
_ x  y _ f  B _ x  y _ 


w  (k  ,k  )  »  - - 

P  x  y  2 

m  w 
P 


-(kx<>  J 


B; 


yi  j 


(11-46) 
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III.  PASS  BAND  CHARACTERISTICS 


In  tht»  analysis  section  we  developed  the  transform  displacement  response 

for  two  infinitely  periodic  structures;  a  plate  supported  by  a  one-dimensional 

array  of  beams,  and  a  plate  supported  by  a  two-dimensional  square  grid  of 

beams.  In  both  cases  the  structure  was  driven  by  a  point  force  located  at  a 

beam  support  location,  and  in  the  grid  case  at  a  grid  intersection  point  as 

well.  One  of  the  purposes  of  this  analysis  was  to  determine  if  such  structures 

exhibited  characteristic  "pass-bands"  within  which  energy  freely  propagated 

from  the  drive  point  throughout  the  structure.  In  a  complementary  sense  the 

determination  of  the  existence  of  non-propagating  bands,  within  which  energy 

was  trapped  in  the  bays  adjacent  to  the  drive  point,  was  also  desired.  The 

purpose  of  this  section  will  be  to  determine  if  such  band  characteristics 

exist  for  the  two  classes  of  structures  investigated.  This  is  accomplished 

by  investigating  the  singularities  in  the  wave  number  spectrum  of  the  plating 

response.  A  given  frequency  will  be  considered  to  be  within  a  pass  band  if, 

at  that  frequency,  singularities  exist  which  fall  on  the  real  (k  ,k  )  axes 

x  y 

thus  representing  undamped  propagating  waves.  Examining  the  limiting 


solution  (k  ,k  <<  k  )  as  represented  by  Eq.  11-35,  we  have 
x  V  t 


w  (k  ,k  )  - 
p  x  y 


I 


I 


ro  u) 
P 


)  v  r/ky  av  ii 
( >  -  f  giv  l(fB)  *  (a)  -  >J ) 


(in-i) 


G(kf) 


cot 


kfL 


+  coth 


¥] 


(IH-2) 


The  limits  of  the  propagation  bands  corresponding  to  this  wavenumber  region 
can  be  determined  from  the  poles  of  Eq.  III-l.  Thus,  examining  the  denominator 
we  see  that  for  a  given  value  of  k^,  we  can  determine  the  zeroes  in  terms  of 
k^,  or  setting 
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4  4 

\  "  2BKB 


B  >  0 


( 1 1 1-3) 


then  we  have 


A 

1/4  ’ 

k  -  ± 

X 

B 

1  1/4 


Y  G(kJ 
m  f 


+  1-B 


(III-4) 


now,  if  the  quantity  within  the  brackets  of  Eq.  III-4  is  positive  then  k^ 
will  be  wholly  real,  implying  a  freely  propagating  wave.  (Alternatively  if 
the  term  within  brackets  is  negative  then  clearly. 


k  -  ±  — .  . . 
x  _1/ 4 


(1+i) 


\nG(V 


+  1-B 


1/4 


( 1 1 1  —  5 ) 


and  the  wave  will  be  damped. 

For  small  values  of  y  (1-B)  the  passband  limits  are  closely  approximated 

in 

by  the  zeroes  of  (G(kf) ]  1,  or 

k,L 


propagating  bands,  nn  <  — n,  n  -  0,1, 2, 3, . . . ,« 


nonpropagating  bands,  n  <  — ~  <  (n+l)n 


0,1, 2, 3, . . . ,® 


For  other  values  of  y  (B-l)  the  plate  wavenumber  limits  of  the  passbands  are 

in 

best  described  graphically.  This  result  is  depicted  in  Figure  III-l  for 

several  representative  values  of  the  beam  to  plate  mass  ratio  parameter,  y  . 

Thus,  the  left  hand  graph  is  a  plot  of  y  (B-l)  vs.  B  for  values  of  y  =  1,5, 

m  ^  in 

and  10.  The  right  hand  graph  is  a  plot  of  [Gikf) ]  vs.  kfL/2.  As  an 

illustrative  example  let  us  assume  that  we  have  a  structure  with  y  =  1. 

m 

Examining  Fig.  III-l  we  see  that  for  ky  wavenumbers  much  less  than  the 


flexural  beam  wavenumber  k_,  the  quantity  y  (B-l)  is  negative  and  approximately 

B  m 

equal  to  -1,  thus  real 
except  those  at  which 


equal  to  -1,  thus  real  (propagating)  values  of  k^  will  result  at  all  frequencies 
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[G(kf)]_1  <  -  1 


Using  the  right  hand  chart  then  we  see  that  this  result  corresponds  to  a 
series  of  narrow  non-propagating  bands,  each  with  a  lower  bound  of 


kfL 


(• 


4n+3^\ 
- )  * 


For  large  values  of 


kfL 


and  values  of 


[y  (B-l) ]  <  -  1 
m  — 

we  can  determine  an  approximate  value  of  the  bandwidth,  from  the  expression 


cot 


NP 


*  *  {*¥)  ■ 


kfL 


r  Ym(B-D 


-  1 


or 


.NP 


As  k  values  nearer  to  k  are  chosen  we  see  that  the  bandwidth  of  the  non- 

y  B 

propagating  bands  increases  until  for. 


Y_(B-1)  >  1 
'  m 

most  of  the  frequencies  correspond  to  non-propagating  bands  with  narrow 
bands  of  propagating  k^  values,  each  with  am  upper  limit  of 


The  bandwidth  in  these  cases  is  given  by 

AP _ 1 _ 

\  ~  1  +  krLy  (B-l) 
t  m 
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The  essential  result  of  this  analysis  is  that  for  any  given  frequency 


there  always  exists  some  set  of  (k  ,k  ) 

x  y 


wavenumbers  which  freely  propagate. 


This  result  has  been  determined  from  an  analysis  of  the  limiting  result  for 


kf  >>  k  ,k  .  Examining  the  other  limiting  results,  as  shown  in  Fig.  II-2, 

•  x  y 

we  see  that  none  of  the  expressions  shown  indicate  passband  characteristics. 


With  a  one-dimensional  structure,  however,  the  results  are  somewhat 
different  in  that  a  well  defined  set  of  pass-bands  does  exist.  The  comparable 
result  from  our  analysis  is  that  provided  by  the  solution  for  the  point  driven 
one-dimensional  array  as  given  by  Eq.  11-44, 


As  th«  grid  problem,  the  limits  of  the  propagation  bands  are  provided  by 
the  poles  of  this  expression,  or  when 


When  the  quantity  in  brackets  is  negative  we  will  have  non-propagating 
wavenumbers  given  by 


The  frequencies  at  which  these  will  occur  are  provided  by  the  expression 

1  \n 

G(kf)  *  2 


and  from  the  result  for  the  grid  we  can  write  with  minor  modification. 
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A.  Propagating  Bands 


nn  < 


2 

«  i 

4n+3 

2  -  \ 

,  4 

,  1m 
2 

»  1 

)  n'  n  =  0, 1, 2, . . .  ,< 


„  V  3, 


and 


Mn+3> 

|  it  +  | 

f  1  ^ 

\ <  V:<( 

140+3^ 

V  4  ) 

kfLy  , 

V1*  2  ) 

r  2 

l  4  / 

B.  Non-Propagating  Bands 

y_ 


1.  —  «  1 


<  (n+l)u  ,  n  =  0,1,2,...,® 


2.  T»1 


^+3^*  kfL  /4n+3\ 


i  + 


k  Ly 
f  m 


2  J 


n  =  0,1,2, ... ,« 


In  comparing  this  result  with  the  previously  cited  work  by  Mead  for  a 
one-dimensional  beam  with  a  periodic  array  of  rotational  restraints,  it  is 
noted  that  the  pass-bands  just  stated  for  the  plate  beam  array  structure 
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compare  with  his  result  for  a  low  rotational  inertia — thus,  an  approximation 

5 

of  a  simply  supported  boundary  condition.  In  this  case  he  shows  that  the 

propagation  bands  are  bounded  by  the  n^  simply  supported  resonance  on  the 

low  side  and  the  nth  clamped-clamped  resonance  on  the  upper  side.  Clearly 

this  is  very  close  to  the  result  just  obtained  in  this  analysis  for  the 

case  y  =0. 
m 


A  further  comparison  with  Mead  results  from  a  comparison  with  the 

solution  for  spring-like  as  opposed  to  mass-like  constraints.  In  the  former 

case  the  pass-bands  flip  with  the  result  that  the  clamped-clamped  resonances 

now  determine  the  lower  boundary  of  the  pass-bands.  An  analogous  behavior 

with  the  grid  supported  plate  is  seen  when  comparing  the  pass  bands  in  either 

low  or  high  (k^/kg)  regions.  Thus  at  low  (k^/kg)  values  the  structure 

represents  an  inertial-like  impedance  to  waves  propagating  in  the  x-direction, 

and  at  high  k  /k  )  values  a  compliance  like  impedance.  This  analogy  becomes 
y  b 

somewhat  more  distinct  if  we  examine  the  behavior  of  the  y  (B-l)  term  in 

'm 


Fig.  III-l.  Thus,  as  k^  values  from  a  small  fraction  of  k&  to  a  large 
multiple  of  k  are  traversed  we  see  a  clear  shift  in  the  sign  of  the  drive 
point  impedance  comparable  in  a  qualitative  sense  to  a  switch  from  an  inertial 
to  a  compliant  reactance.  It  is  this  periodic  shift  which  directly  results 
in  the  variance  between  the  one  dimensional  and  two  dimensional  pass-band 
results.  The  relevance  of  the  above  results  to  structure  borne  noise  problems 
is  as  follows.  For  the  one  dimensional  case,  the  existence  of  finite  band 
width  pass  bands  and  stop  bands  implies  that  such  a  structure  can  be  used 
as  a  natural  filter  to  attenuate  structure  borne  noise  propagation  from  a 
localized  excitation.  In  other  words,  the  one  dimensional  structure  may,  in 
theory,  be  tuned  so  that  particularly  offensive  portions  of  the  frequency 
spectrum  of  the  excitation  correspond  to  stop  bands  of  the  structure.  However, 
the  two  dimensional,  or  grid,  situation  is  different.  Here,  there  are  no 
stop  bands  if  one  assumes  that  the  vibration  source  exhibits  a  broad  wave- 
number  spectrum  such  as  is  produced  by  a  localized  vibration  source.  In 
other  words  there  is  always  some  portion  of  the  wavenumber  spectrum  of  the 
excitation  that  corresponds  to  a  pass-band  and  thus  short  circuits  the  stop 
band  phenomenon.  Therefore,  in  this  case,  in  addition  to  tuning  the  frequency 
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response  of  the  structure  to  the  frequency  spectrum  of  the  excitation  one 
must  also  tune  the  wavenumber  response  of  the  structure  to  the  wavenumber 
spectrum  of  the  source.  Since  the  wavenumber  spectrum  of  a  source  is 
determined  largely  by  the  geometrical  configuration  of  its  foundation 
including  the  structural  details  of  the  footings,  any  attempt  to  tune  a 
grid  structure  to  take  advantage  of  the  existence  of  stop  bands  would 
have  to  consider  such  design  details. 


-22- 


io  *5  j^reerE-T) 


CAMBRIDGE  ACOUSTICAL  ASSOCIATES.  INC 


IV.  RADIATED  PRESSURE 

An  interesting  and  important  result  directly  obtained  from  the  analysis 
of  Section  II  is  the  radiated  sound  pressure  measured  at  a  point,  (R,0,$),  0 
being  measured  from  the  normal  to  the  grid  or  array  supported  plate  and 
passing  through  the  drive  point  location.  This  result  can  be  written  as,6 


p  ( R,  0 , 4>)  =  -  poo  -  w(ksin0cos<J>,  ksin0sin$) 


(IV-1) 


Thus,  if  fluid  loading  is  neglected  the  far  field  radiated  pressure  can  be 
developed  in  a  closed-form  using  Eqs.  11-27  and  11-32.  In  the  balance  of 
this  section  we  will  examine  the  effect  of  various  structural  configurations 
on  the  "on-axis"  radiated  pressure. 


p(R,0,0) 


w(0,0) 


where  w(0,0)  is  given  in  Eq.  11-36.  We  will  compare  the  resultant  radiated 
pressure  from  four  different  structural  configurations.  In  making  these 
comparisons  we  have  used  a  reference  pressure,  pQ,  equivalent  to  the  high 
frequency  limit  for  the  point  driven  infinite  plate. 


pF  ikR 
o  e 


(IV-2) 


A.  Point  Driven  Infinite  Plate 


This  well  known  result  is  given  by, 


P  (R,0,0) 


(IV- 3) 


a)  -►  0 


(IV-4) 
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B.  Point  Driven  Array  Supported  Plate 

Two  results  are  available  here,  one  for  the  drivepoint  at  a  midspan 
location,7  and  the  one  provided  by  Eq.  11-45  for  a  drivepoint  at  a  beam 
support  location. 


1.  Mid  Span  Drive  Point 


pam(r'°'0) 

P  (R,0,0) 

oo 

K- 

sinkfL/2  '  sinhkfL/2  | 

Po 

Po 

)  - 

2  kfL  kfL  ( 

k  a  +  cot  ~r +  coth  — ' 

\  ( IV-6) 


[t]  '  “  *  ° 


(IV-7) 


(¥)—(¥) ' 


(IV-8) 


2.  Beam  Drive  Point 

From  Eq.  11-45  we  can  write  directly, 


Pab(r,°,0)  p~( 


r,o,o) |  r 

To  7  V  kfL  S"!  (IV"9) 

L  1  +  —  —  (cot  —  +  coth  —  jj 


k]  • 


io  -►  0 


(IV- 10) 


v”  kfL  r  .  k«L  ,i 

7  T  [=ot  -3-  +  ij  ■ 


(0  •+  * 


(IV-12) 
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C.  Point  Driven  Grid  Supported  Plate 

As  stated  in  the  analysis  section  the  radiated  pressure  is  developed 
from  Eq.  11-36.  It  should  benoted  that  in  contrast  to  the  results  presented 
in  the  analysis  section,  where  the  limiting  solutions  were  developed  with 
fluid  loading  neglected,  we  have  here  retained  fluid  loading  to  the  extent 
that  it  is  represented  in  the  low  frequency  limit  of  the  solution.  Thus, 
from  Eq.  II-8, 


yp(0,0) 


P  (R,0,0) 


T  2 
L<o  m 


(IV-13) 


Using  the  transform  structural  displacement  from  Eq.  11-36  we  have  for  the 
radiated  pressure. 


(IV-14) 


( IV-15) 


V 


(cot  -f-  +  l) 


(IV-16) 


Summaries  of  these  results  are  shown  in  Figs.  IV-1  and  IV-2.  The  dimensions 

are  typical  of  airframe  structures.  A  structural  damping  factor  of  np  •  .1  was 

also  used.  Several  points  are  of  interest  to  this  analysis.  At  low 

frequencies  the  point  driven  structures  differ  from  the  infinite  plate  result 

by  a  factor  determined  solely  by  the  ratio  of  beam  to  plate  inertia,  v  .  The 

in 

actual  amounts  are  shown  in  each  figure. 

At  high  frequencies  the  point  driven  structures  fall  off  as  20  log  (f) 
for  the  beam  drive  point  cases,  and  approach  the  infinite  plate  result,  which 
is  frequency  independent,  for  the  midspan  drive  point  case.  This  is  due  to 
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the  fact  that  the  drive  point  impedance  of  a  point  driven  infinite  beam 
is  proportional  to  frequency  while  that  of  the  infinite  plates  is  frequency 
independent. 


A  further  point  of  interest  is  provided  by  a  comparison  with  the  pass 
band  characteristics  determined  in  Section  III.  The  peaks  and  valleys  in 
Figs.  IV-1  and  IV-2  are  related  to  the  interframe  plating  resonances  and 
antiresonances  respectively.  These  plating  resonances  imply  large  reactions 
at  the  supports  adjacent  to  the  drive  point  which  serve  to  reflect,  rather 
than  transmit,  incoming  waves.  Thus  the  bands  corresponding  to  efficient 
levels  of  radiation  are  very  close  to  the  non-propagating  bands  of  the 
structural  analysis. 
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